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A fast and efficient numerical-analytical approach is proposed for modeling complex behaviour in the BBGKY 
hierarchy of kinetic equations. We construct the multiscale representation for hierarchy of reduced distribution 
functions in the variational approach and multiresolution decomposition in polynomial tensor algebras of high- 
localized states. Numerical modeling shows the creation of various internal structures from localized modes, which 
are related to localized or chaotic type of behaviour and the corresponding patterns (waveletons) formation. The 
localized pattern is a model for energy confinement state (fusion) in plasma. 



1. INTRODUCTION 

Kinetic theory is an important part of general 
statistical physics related to phenomena which 
cannot be understood on the thermodynamic 
or fluid models level [1]. In these two papers 
we consider the applications of a new numeri- 
cal/analytical technique based on wavelet anal- 
ysis approach for calculations related to the de- 
scription of complex (non-equilibrium) behaviour 
of the corresponding classical and quantum en- 
sembles. The classical ensembles in this part 
are considered in the framework of the general 
BBGKY hierarchy and the quantum ones in part 
2 in the Wigner-Weyl approach [1]. We restrict 
ourselves to the rational/polynomial type of non- 
linearities (with respect to the set of all dynami- 
cal variables) that allows to use our results from 
[2], which are based on the so called multireso- 
lution framework [3] and the variational formula- 
tion of initial nonlinear (pseudodifferential) prob- 
lems. Wavelet analysis is a set of mathematical 
methods which give a possibility to work with 
well-localized bases in functional spaces and pro- 
vide the maximum sparse forms for the general 
type of operators (differential, integral, pseudod- 
ifferential) in such bases. It provides the best pos- 
sible rates of convergence and minimal complexity 
of algorithms inside and, as a result, saves CPU 



time and HDD space [3]. Our main goals are 
an attempt of classification and construction of a 
possible zoo of nontrivial (meta) stable states: (a) 
high-localized (nonlinear) eigenmodes, (b) com- 
plex (chaotic-like or entangled) patterns, (c) lo- 
calized (stable) patterns (waveletons). In case (c) 
an energy is distributed during some time (suffi- 
ciently large) between only a few localized modes 
(from point (a)). We believe, it is a good image 
for plasma in a fusion state (energy confinement). 
Our construction of cut-off of the infinite system 
of equations is based on some criterion of conver- 
gence of the full solution. This criterion is based 
on a natural norm in the proper functional space, 
which takes into account (non-pcrturbatively) the 
underlying multiscale structure of complex statis- 
tical dynamics. In Sec. 2 the kinetic BBGKY 
hierarchy is formulated. In Sec. 3 we present 
the explicit analytical construction of solutions 
of the hierarchy, which is based on tensor alge- 
bra extensions of bases generated by the hidden 
multiresolution structure and proper variational 
formulation leading to an algebraic parametriza- 
tion of the solutions. So, our approach resem- 
bles Bogolyubov's and related approaches but we 
don't use any perturbation technique (like virial 
expansion) or linearization procedures. Numeri- 
cal modeling as in general case as in particular 



cases of the Vlasov-like equations shows the cre- 
ation of various internal structures from locaUzed 
bases modes, which demonstrate the possibhty of 
existence of (metastable) pattern formation. 

2. BBGKY HIERARCHY 

Let M be the phase space of an ensemble 
of N particles (dimM = 6-/V) with coordinates 

Xi = {qt,Pi), i = l,---,^, qi = {q},qf,qf) G 

^^ Pt^iphphpl) <^R^, q = {qu---,qN) e 

R^^ . Individual and collective measures are: 
fit = dxt = dq^dpt, fi = rii^iMi- Our con- 
structions can be applied to the following general 
Hamiltonians: 

N 2 

^" = E(|^ + ^«(9))+ E U,,{q„qj) 
i=l l<i<j<N 

where the potentials Ui{q) — Uiiqi, . . . , qN) and 
Uij{qi,qj) are restricted to rational functions of 
the coordinates. Let Ls and L^ be the Liouvil- 
lean operators (vector fields) 






Pl^ 

m dqj dq dpj 



dUj d 



E L^^ (1) 

l<i<j<s 



*^ dqi dp.. 



9C/y d 



dq-j dpj 

and Fn{xi, . . . ,XN',t) be the hierarchy of A^- 
particle distribution function, satisfying the stan- 
dard BBGKY-hierarchy (V is the volume) [1]: 



dt 



+ LsFs 



1 



d/^s+l ^2^ii,s + l-Fs+l (2) 



In most cases, one is interested in a representation 
of the form Fk{xi, . . . ,Xk;t) = lli=iFi{xi;t) + 
Gk{xi, . . . jXk'jt) where Gk are correlators. Ad- 
ditional reductions often lead to simplifications, 
the simplest one, Gk — 0, corresponding to the 
Vlasov approximation. Such physically moti- 
vated ansatzcs for Fk formally replace the linear 
(in Fk) and pseudodifferential (in general case) 
infinite system (2) by a finite-dimensional but 
nonlinear system with polynomial nonlinearities 
(more exactly, multilinearities [3]). Our key point 
in the following consideration is the proper non- 
perturbative generalization of the perturbativc 
multiscale approach of Bogolyubov. 



3. MULTISCALE ANALYSIS 

The infinite hierarchy of distribution functions 
satisfying system (2) in the thcrmodynamical 
limit is: 

F = {Fo, Fi(xi; i), . . . , Fn{xi, ...,XN]t),.. .} (3) 

where Fp{xi, . . . ,Xp]t) e HP, H" = R, HP ^ 
L^{R^P) (or any different proper functional spa- 
ce), F e H°° = H^ ®H^ ® ...®HP ® ... with 
the natural Fock space like norm (guaranteeing 
the positivity of the full measure) : 



{F,F)=FS+Y,jFi{x^ 



;^)n 



A^^ (4) 



£=1 



First of all we consider F — F{t) as a function 
of time only, F E L^ (R) , via multiresolution de- 
composition which naturally and efficiently intro- 
duces the infinite sequence of the underlying hid- 
den scales [3] . Because the aflfine group of transla- 
tions and dilations generates multiresolution ap- 
proach, this method resembles the action of a mi- 
croscope. We have the contribution to the final 
result from each scale of resolution from the whole 
infinite scale of spaces. We consider a multireso- 
lution decomposition of L'^{R) [3] (of course, we 
may consider any different and proper for some 
particular case functional space) which is a se- 
quence of increasing closed subspaces Vj G L^ (R) 
(subspaces for modes with fixed dilation value): 



.V-2 C V-i C Fo C T/i C ^2 C 



(5) 



The closed subspace Vj{j <E Z) corresponds to 
the level j of resolution, or to the scale j and 
satisfies the following properties: let Wj be the 
orthonormal complement of Vj with respect to 
Vj+i- Vj+i = ^ ^j- Then we have the follow- 
ing decomposition: 



{Fit)} 



w,=v,^w, 



(6) 



- oo <j < oo 



j=o 



in case when Vc is the coarsest scale of resolution. 
The subgroup of translations generates a basis 
for the fixed scale number: spa.nf.^^{2^^^'^{2H — 
k)} — Wj. The whole basis is generated by action 



of the full afFine group: 

spanfcez,,ez{2-'/2vl/(2^t_fc)}= (7) 

spanfc,,e2{*j,fc} = {F{t)} 

Let the sequence {Vj},Vj C L'^{R) corresp- 
ond to multiresolution analysis on the time ax- 
is, {V.^'} correspond to multiresolution anal- 
ysis for coordinate Xi, then V"'^^ = V^^ ® 
. . . ® V.^" (g) Vj corresponds to the multires- 
olution analysis for the n-particle distribution 
function i^„(xi, . . . , a;„; i). E.g., for n = 2: 
Vo = {/ : fixi,X2) = Efci,fc2"fei^fe2'A^(2;i - 
fci, 2:2-^2), afei,fc2 e •^^(2'2)}, where 0^(xi, 2:2) = 
(I)^{xi)(f)'^{x2) = (j)^ ® (/)2(xi,a;2), and ^'(xi) = 
(j){xi) form a multiresolution basis corresponding 
to {Vp}. If {(t>^{xi - t)], I ^ Z form an or- 
thonormal set, then (jP'(xi — ki,X2 — fe) form 
an orthonormal basis for Vq. So, the action of 
the affine group generates multiresolution repre- 
sentation of L^(i?^). After introducing the de- 
tail spaces W^, we have, e.g. V^ = Vq ® Wq. 
Then the 3-component basis for W§ is gener- 
ated by the translations of three functions: ^^ = 
0i(a;i) ® *2(a;2)^^2 ^ <f i(xi) » <?!)2(x2), *§ = 
^^(xi) (g) ^^(x2). Also, we may use the rectan- 
gle lattice of scales and one-dimensional wavelet 
decomposition: 

f{xi,X2)= Y^ (/,*j,£®*j,fe)*j,£«)*j,fc(a;i,a;2) 

i,£:j.k 

where the basis functions '^i^i "^j.k depend 
on two scales 2^' and 2^^. We obtain our 
multiscale/multiresolution representations (for- 
mulae (11) below) via the variational wavelet ap- 
proach for the following formal representation of 
the BBGKY system (9) (or its finite-dimensional 
nonlinear approximation for the n-particle dis- 
tribution functions) with the corresponding ob- 
vious constraints on the distribution functions. 
Let L be an arbitrary (non) linear differential- 
/integral operator with matrix dimension d (fi- 
nite or infinite), which acts on some set of func- 
tions from L2(Vj®"). -qj ^ '^{t,Xi,X2, ■ ■ ■) = 

(^*i(i,xi,X2,... ),..., *''(i,xi,a;2,...)), x, e 

Q. C R®, n is the number of particles: 



Q = Qdo,di,d2,...{tTXi,X2, . • . , 

d/dt,d/dxi,d/dx2,..., / ^fc) 



da,di,d2,... 

/ J Q_ioiii2...V'iXi,X2, ■ ■ ■) 

io,n,J2,---=i 

\dtJ \dxiJ \dx2J J 

Let us consider now the N mode approximation 
for the solution as the following ansatz: 



N 



'ii^{t,xi,X2,...) = ^ atotit2... 



(9) 



io,ilA2,--- = l 



Aig ®Bi^ (X)Ci2 ...{t,xi,x2,...) 



We shall determine the expansion coefficients 
from the following conditions (different related 
variational approaches are considered in [2]): 






(10) 



(L*^' )Ak„ {t)Bk, {xi)Ck2 (a;2)dMxida;2 . . 



Thus, we have exactly diV" algebraical equations 



for dN"^ unknowns a^, 



So, the solution is 



parametrized by the solutions of two sets of re- 
duced algebraical problems, one is linear or non- 
linear (depending on the structure of the operator 
L) and the rest are linear problems related to the 
computation of the coefficients of the algebraic 
equations (10). which can be found by using the 
compactly supported wavelet basis functions for 
the expansions (9). As a result the solution of the 
equations (2) has the following multiscale decom- 
position via nonlinear high-localized eigenmodes 

F{t,xi,x2,...)= Yl a,jW(g)V^{t,Xi,...) 



VHt) = V^-'^^it) + J2 V/iu^it), ui ^ 2' (11) 

1>N 

U (Xs) = Uj^ {xs)+ 2_^ U^{k^Xs), k^ ^ 2 



m>M 



L^ = L{Q,t,x,)^{t,x,)^0, 



(8) 



which corresponds to the full multiresolution ex- 
pansion in all underlying time/space scales. The 



formulae (11) give the expansion into a slow and 
fast oscillating parts. So, we may move from the 
coarse scales of resolution to the finest ones for ob- 
taining more detailed information about the dy- 
namical process. In this way one obtains con- 
tributions to the full solution from each scale 
of resolution or each time/space scale or from 
each nonlinear eigenmode. It should be noted 
that such representations give the best possi- 
ble localization properties in the corresponding 
(phase)space/time coordinates. Formulae (11) 
do not use perturbation techniques or lineariza- 
tion procedures. Numerical calculations are 







Figure 3. Localized waveleton pattern. 




Figure 1. Localized mode contribution to distri- 
bution function. 



based on compactly supported wavelets and re- 
lated wavelet families [3] and on evaluation of 
the accuracy on the level N of the correspond- 
ing cut-off of the full system (2) regarding norm 
(4): ||i^^+i - F^W < e. We believe that the 
appearance of nontrivial localized patterns (a)- 
(c) demonstrated on Fig.l-Fig.3 constructed by 
these methods is a general effect which is also 
present in the full BBGKY hierarchy, due to 
its complicated intrinsic multiscale dynamics and 
it depends on neither the cut-off level nor the 
phenomenological-like hypothesis on correlators. 
So, representations like (11) and the prediction of 
the existence of the (asymptotically) stable local- 
ized patterns/states (energy confinement states) 
in BBGKY-like systems are the main results of 
this paper. 
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